Excitons and cavity polaritons for ultracold atoms in an optical lattice 
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We study the resonant electronic excitation dynamics for ultracold atoms trapped in a deep optical 
lattice prepared in a Mott insulator state. Excitons in these artificial crystals are similar to Frenkel 
excitons in Noble atom or molecular crystals. They appear when the atomic excited state line 
width is smaller than the exciton band width generated by dipole-dipole coupling. When the atoms 
are placed within a cavity the electronic excitations and the quantized cavity mode get coupled. 
In the collective strong coupling regime excitations form two branches of cavity polaritons with 
Rabi splitting larger than the atomic and the cavity line width. To demonstrate their properties 
we calculate the transmission, reflection, and absorption spectra for an incident weak probe field, 
which show resonances at the polariton frequencies. 

PACS numbers: 42.50.-p, 71.36.+C, 71.35.Lk 



I. INTRODUCTION 

The quantum phase transition from the superfluid to 
the Mott insulator phase of cold dilute gas of boson atoms 
in an optical lattice was widely studied both theoretically 
[J, 0] , and experimentally [1] . The optical lattice poten- 
tial can be formed of counter propagating laser beams, 
where their interference forms a standing wave with pe- 
riod of A/2, half wave length [4(. The geometry and the 
depth of the optical lattice potential are easily controlled 
by the laser field. The cold atoms loaded on the optical 
lattice realize the Bose-Hubbard model [5|, which pre- 
dicts the phase transition from the superfluid into the 
Mott insulator phase, by increasing the optical potential 
depth. The atoms are confined in an array of microscopic 
trapping potentials due to the optical dipole force, and 
by changing the relevant parameters, the optical lattice 
can be filled with few atoms down to one atom per site 
0, Hi ■ Optical lattices allow us to go beyond the weakly 
interacting regime of Bose gas into strongly correlated 
systems of condensed matter physics [||. 

The cold atoms in an optical lattice, in the Mott in- 
sulator phase with one atom per site, can be considered 
as an artificial crystal. Such artificial crystals are simi- 
lar to Noble atom and molecular crystals, as each atom 
retains his identity, and negligible overlaps exist between 
the atom electronic wave functions at different sites. It 
will be of big importance to check the possibility of the 
appearance of solid state effects and processes in such 
artificial crystals. For example in molecular crystals, an 
electronic excitation localized on a molecule can transfer 
among the crystal molecules due to electrostatic inter- 
actions, i.e. dipole-dipole interactions. Namely, an ex- 
cited molecule will decay to the ground state and other 
molecule will be excited, without a charge transfer be- 
tween the molecules. In using the crystal symmetry, 
such energy transfer can be represented in the quasi- 
momentum space by a wave that propagates in the crys- 
tal. We obtain quasi-particles with an effective mass and 
a finite life time, which are called excitons If the 

molecular crystal is located within a cavity, the excitons 



and the cavity photons are coupled. In the strong cou- 
pling regime, where the exciton and photon line-widths 
are smaller than the coupling strength, the excitons and 
photons are coherently mixed to form new system quasi- 
particles which are called cavity-polaritons [g, Q . We get 
two polariton branches which are separated by the Rabi 
splitting energy. Coherent interactions of laser pulses in a 
resonant optically dense medium of two-level atoms was 
presented in fioj ]. 

In the present paper we investigate such excitons and 
polaritons. We consider two- level atoms within a cav- 
ity, with transition energies close to resonance with a 
cavity mode. The optical lattice potentials are chosen 
to give light shifts in-phase for the ground and excited 
states, with minimums at the same positions. In mul- 
tilevel atoms and for a specific laser frequency one can 
achieve an excited level with light shifts in-phase with 
the ground state ones, and there even exists a magic fre- 
quency where we get equal shifts as in the ground state 
p, [ll[ . In all the discussions the trapping fields are con- 
sidered classical and far from resonance with any of the 
atom electronic transitions. Furthermore, we limit the 
treatment to the Mott insulator phase with one atom 
per site [Hj]. We examine the possibility of the appear- 
ance of excitons in such a system, by comparing between 
the energy transfer coupling parameter induced by the 
dipole-dipole resonance interaction, and the excited state 
line width. In a cavity the view is different, now the ex- 
citons are strongly coupled to the cavity photons by the 
dipole interaction. Hence, coherent polariton states can 
be observed which are superpositions of excitons and cav- 
ity photons. In paper [13|], the cavity photons are also 
treated quantum mechanically, but their energies are far 
to resonance with any of the atom transitions, and they 
are used to produce the optical lattice. Moreover, we get 
the linear optical spectra of the cavity for an incident 
external field. 

The paper is organized as follows. In section 2 we 
discuss excited and ground state cold atoms loaded on 
an optical lattice. The excitons in such a system are 
introduced in section 3. Section 4 describes polaritons in 
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a system of optical lattice cold atoms within a cavity. The 
transmission, reflection, and absorption spectra are given 
in section 5. The conclusions are presented in section 6. 
In the appendix the linear optical spectra are derived 
using the input-output formalism. 



II. EXCITED AND GROUND STATE 
COLD-ATOMS IN AN OPTICAL-LATTICE 

We consider neutral cold two-level atoms where the 
cavity mode is close to resonance with only a single ex- 
cited state. The cold atoms are trapped on 2D standing 
laser waves far to resonance with any of the atomic ex- 
citations, and result in shifts in the ground and excited 
states. Such shifts correspond to an optical lattice poten- 
tial V^(x) for ground state atoms, and V e (x) for excited 
state atoms, where x is the atom location. The potentials 
are periodic with a periodicity of, A/2, half laser wave 
length. The potentials are in general out of phase, and 
for two-level atoms they have ir/2 phase shift, where the 
position of the minimum of one potential appears at the 
maximum of the other. Even though, we assume the two 
optical lattice potentials to be in-phase, with minimums 
located at the same position. Such assumption is possi- 
ble as for multi-level atoms, i.e. Alkali or Alkaline atoms, 
one get a number of optical lattices, one for each level. 
For a specific laser magic frequency and polarization ex- 
ists an excited state optical lattice potential is equal to 
the ground state one [y, [ll[ . The two dimensional opti- 
cal lattice potentials are V s (x) = Yli^s cos 2 (kxi), with 
(xi — x.y), and (s = g,e), where k — 2ir/\, and V s is 
the s potential depth, see figure (1). Furthermore, the 
ground and excited state atoms are located in the first 
Bloch band of the ground and excited optical lattice po- 
tentials, which is the case for deep optical lattices. 




Ground state optical lattice 



FIG. 1: Schematic plot of the ground and excited state optical 
potentials. 

The ground and excited state atoms are considered as 
two kinds of bosons, which are represented by the Bosc- 
Hubbard model of two components of bosons [l2>], that 
yields a rich phase diagram with two main phases, the 
Mott insulator and the superfluid phases. In the present 
paper we consider only the Mott insulator phase with one 
atom per site, which was realized experimentally Q. 



The internal two-state Hamiltonian of a free atom is 
H A = tkj' g \g')(g'\ + tkj' e \e')(e'\, where \g') and |e') are the 
ground and excited states with energies huj' g and huj' e , 
respectively. Now the atoms are loaded on an optical 
lattice, and therefore the atom internal states become a 
function of the atom position on the optical lattice. In 
the Mott insulator phase with one atom per site, and for 
the first Bloch band, the Hamiltonian of the atoms reads 

H a = ^ i huj a\^)(9i\ + fi^e\e l }(e l \} . (1) 

i 

The summation is over the optical lattice sites. The opti- 
cal lattice in this case results, first in energy shifts of the 
atom internal states, where we have now the energies hoj g 
and ftw e in place of the free atom ones, and the atomic 
energies are the same at each site. At the laser magic 
frequency the atom frequency transition uj e — ujf equals 
to the free atom transition u>' e — lu'^. Second, the optical 
lattice limits the atom positions to be confined on a lat- 
tice with lattice constant a = A/2, which is indicated by 
the index i in the above Hamiltonian. 

Each atom retains his identity. Namely, as the average 
distance between the atoms is of the order of the lattice 
constant a, no overlaps of the electronic wave functions 
exist between atoms at different sites. This fact makes 
such artificial atomic lattices similar to Nobel atom or 
Molecular crystals. Here the crystal forms of external 
laser fields which can be controlled, e.g. it is possible to 
control the lattice constant and the number of atoms per 
site. Even though, the electronic excitations can transfer 
between atoms at different sites due to electrostatic in- 
teractions between the atoms, i.e. dipole-dipole interac- 
tions, as an excited atom has a transition electric dipole. 
An excited atom at site i jumps to the ground sate and 
another atom at site j is excited, no charge transfer exists 
between the atoms of such interactions, see figure (2) for 
ID case. The resonance energy transfer can be induced 
by dipole-dipole interactions, or higher order multi-pole 
interactions, with the energy transfer parameter Jy . The 
energy transfer Hamiltonian between atoms is given by 

H T = Y,^H S\S h (2) 
id 

where Sj = \ei){gi\ and Si = |<ft)(ej|, are creation and 
annihilation operators of an excitation in the atom at 
site i. 

The atomic Hamiltonian, in terms of internal excita- 
tion creation and annihilation operators, is given by 

H a = J2 Si Si + hJ v S i S h (3) 

i id 

where cu a — lu € — uj g . The Hamiltonian includes only 
terms that conserve the excitation number. 
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FIG. 2: Energy transfer in ID atoms chain. 

III. EXCITONS IN OPTICAL LATTICE COLD 
ATOMS 

The electronic excitations can transfer among the 
atoms in the optical lattice with energy transfer param- 
eter Jjj, between two atoms at sites % and j. The cou- 
pling parameter is a function of the distance between the 
atoms, Jij = J(|rii — rij|), where n, is the atom i location 
in the lattice. The excitation can be at any atom with 
the same probability, and in using the lattice symme- 
try the excitation is represented in the quasi-momentum 
space by a wave propagates in the lattice with wave vec- 
tor k. Such a quasi-particle is called Frenkel-exciton in 
molecular crystals, e.g. in Nobel atom or organic crystals 
0, Q • Here we study Frenkel like-excitons in optical lat- 
tice cold atoms, and we give the condition for obtaining 
such quasi-particles in such a system. 

The operators Si are of two-level atom which are as for 
a system of spin-half lattice. They obey Bose commuta- 
tion relation between different sites, namely [Si, Sj] = 0, 
and Fermi anti-commutation relation on the same site, 
namely {Si,Sj} = 1, which forbids two excitations from 
being localized on the same atom. At low excitation 
density we neglect the possibility of getting two exci- 
tations on the same atom, and we can assume the ex- 
citations to behave as Bosons 14]. The bosonization 
procedure is done by making the replacement Si — > Bi, 
where [Bi,Bj] — 5ij. Namely we neglect saturation ef- 
fects which give rise to nonlinear processes. We treat 
the single excitation case, which is given by the bosonic 
Hamiltonian 

i i,i 

The Hamiltonian can be easily diagonalize by the trans- 
formation into the quasi-momentum space, in using 

where N is the number of lattice sites. The exciton 
Hamiltonian reads 

fl- a = 5^fiw a (k)BiB k , (6) 

k 



with the exciton dispersion 

a, a (k)=a, a + ^;j(L) e-' ;k ' L , (7) 

L 

where we defined the distance between atoms by L = 
rii - rij. 

Our system is 2D square with a cubic symmetry. 
The in-plane wave vector is defined by k = (k x ,k y ) = 
(n X: n y ) x 27r/v / 5, where the system area is S = Na 2 , 
and we have (n x ,n y — 0, ±1, ••• ,±y/~N/2). If we as- 
sume energy transfer only between the first two nearest 
neighbor sites with coupling parameters J(a) = —J\ and 
J(v2a) = — Ji, which are negative and give attractive 
interactions, hence the exciton dispersion reads 

w a (k) = uj a — 2J\ (cos k x a + cos k y a) 

— 4J 2 cos(V2ak x ) cos(V2ak y ). (8) 

The energy transfer in a symmetric lattice results in an 
energy band, in place of a discrete level for independent 
atoms. In considering only the first nearest neighbor in- 
teraction, the band width is 47iJi. In fact the dipolc- 
dipole interaction is of long-range interaction, and hence 
one needs to go beyond the nearest neighbor interactions. 
But here as the lattice constant is of the order of the 
atom transition wave length, we limited the discussion 
for nearest neighbor interactions. 

In the long wave length limit, that is ka <C 1 with k = 
|k|, which is the regime of importance for the excitation- 
photon coupling as we will see in the next section, we get 
the parabolic dispersion 

hk 2 

uj a (k) = cu a - 4( Ji + J 2 ) + , (9) 

Zm eff 

where we defined the exciton effective mass by m e // = 
h/[2a 2 (Ji + 4J 2 )]. Note that at zero wave vector, k — 0, 
a significant atomic frequency shift of 4( J\ + J2), relative 
to the free atom frequency, is obtained and which can be 
easily observed. Also such a shift has influence on optical 
lattice clocks (lEj . 

The exciton coherent states defined in Eq.© compete 
with the excited state life time r a . In order to observe 
exciton effects in optical lattice cold atoms, the exciton 
life time need to be longer than the energy transfer time 
1/ Ji between nearest neighbor sites; or the excitation line 
width, hj a — h/r a , need to be smaller than the exciton 
band width of the nearest neighbor interaction, namely 
^-7a < 4hJ\. For example, in free 85 Kb for the transi- 
tion 5 2 Si/2 — 5 2 P3/2 we have hu> a — 1.56 eV, with line 
width of fvj a = 2.5 x 10~ 8 eV; and for free 23 Na for 
the transition 3 2 Si/ 2 — 3 2 P 3 / 2 we have hu a = 2.1 eV, 
with h-f a = 4x 10~ 8 eV. The transfer energies hJ\ and 
TiJi are calculated from the dipole-dipole interaction be- 
tween two atom transition dipoles pi 1 and fx 2 which are 
separated by the distance vector R — RR, which is given 
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by 

1 2 

hJ{R) = V { (<% - 3^41 (cos Zi? + Zi? sin Zi?) 

* — ' 47reo-ft J IV / 

hi 

- Uij-RiRj) l 2 R 2 coslRj, (10) 

here — x,y, z), with oj a — cl. We assume R = RR Z , 
and for the case of two equal dipoles in the z direction, 
where fi 1 = fl 2 = [iR zl we get the coupling 

u 2 

hJ(R) = - — ^— (cos IR + IR sin IR) . (11) 

The distance between two atoms, for the nearest neighbor 
coupling, equals the lattice constant, that is R = a. For 
atoms with hw a — 2 eV, dipole of /i = 2 eA, and optical 
lattice of a = 1000 A, we get HJ\ = 10 -7 eV, which is 
one order of magnitude larger than the free atom line 
width, and excitons can be observed in this case. For 
the next nearest neighbor interaction we get hJ^ = 3 x 
10~ 8 eV, which is of the order of the free atom line width. 
The higher nearest neighbor transfer interaction terms 
are much smaller than the free atom line width, and they 
are not relevant for the exciton formation, and hence can 
be neglected here. 

In the next section we show how in using a cavity one 
can observe better coherent effects involving cavity pho- 
tons, where the excitation transfer is mediated by the 
cavity photons, in the strong coupling regime, where the 
excitation-photon coupling is larger than the atomic and 
cavity-photon line widths. 



IV. CAVITY POLARITONS IN A COLD ATOM 
CRYSTAL 

Let us now consider such a quantum gas of ultracold 
atoms enclosed in the middle of optical resonator built 
of two parallel mirrors, which we consider as perfect for 
the moment (see figure (3)). Note that in contrast to 
most cavity QED configurations, where one considers sin- 
gle modes of spherical mirror resonators with Gaussian 
transverse geometry, our cavity is adapted to the trans- 
lational invariance of the lattice supporting plane wave 
field modes. 

As the electromagnetic field is not confined in the cav- 
ity plane with in-plane wave vector q, the modes have dis- 
crete wave vector k z = rrm/L, where m = 1, 2, • • ■ only 
in the perpendicular z direction. As the cavity distance 
in z direction is small we can restrict our considerations 
to only the m-mode close to resonance with the atomic 
electronic excitation. As the optical lattice is located in 
the middle between the cavity mirrors the active mode 
is taken to be one of the odd modes to ensure strong 
coupling. For each in-plane wave vector exist two possi- 
ble polarizations, TE and TM modes. Here we consider 
only isotropic materials so that we can concentrate in a 




FIG. 3: Optical lattice cold atoms within a cavity. 



fixed cavity photon polarization adapted to the atomic 
excitations. The cavity Hamiltonian thus is given by [l7j 

H c =^ huJ c(q) OqOq, (12) 

q 

where a q , a q are the creation and annihilation boson 
operators of a cavity photon with in-plane wave vector q. 
The corresponding cavity photon dispersion thus reads 




where L is the distance between the cavity mirrors, and 
e « 1 is the cavity medium dielectric constant. Note the 
following discussions also hold for the case of closed cavi- 
ties with discrete and completely confined modes. Single 
mode operation corresponds to limiting ourselves to the 
case of zero in-plane wave vector, k — 0, and with a 
fixed detuning between the atom transition and the cav- 
ity mode. 

The atom electronic transitions and thus the excitons 
are coupled to the cavity modes by the dipole interaction 

V = —jx ■ E, where fi = + &i) i s trie dipole 

operator and E is the cavity field operator. Assuming 
equal dipole moments for all the atoms and the Mott 
insulator phase with one atom per site, the interaction 
Hamiltonian in the rotating wave approximation reads: 

v = E W k * s l ak + /w 4$} , (14) 

ki 

where the coupling parameter is 

^ = (/ ^ k)eIk ' n '' (15) 

S is the mirror area, and is the photon polarization 
unit vector. In using the previous bosonization proce- 
dure, and in applying the transformation (|5|), we get 

v = E h (/* s k a k + k « k ^k) , (i6) 

k 
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with 



excitonic, that is \X k | 2 



Mk 



ltkJ c (k)Nii 2 
2LSe 



(17) 



Here we assumed an isotropic atomic ground state so 
that for each direction of k we have the same dipole 
moment /i. The atom position dependence disappears 
due to the lattice translational symmetry, and the lattice 
symmetry ensures that only coupling between excitons 
and photons with the same in-plane wave vector occurs, 
i.e. we have quasi-momentum conservation. The coupled 
exciton-photon Hamiltonian is written as 

H = ^2h{u a (k) B k B k + w c (fc) 4a k 

k 

+ fk 5 k a k + ft a k 5 k } . (18) 

Due to exciton-photon coupling the exciton and the 
photon coherently mix and it makes sense to transform 
to a diagonal basis with respect to this coupling, namely 
the polariton basis [8, 9]. The maximal coupling appears 
at small wave vectors as the exciton band width is nar- 
row relative to the photon band and the exciton effective 
mass thus is much larger than the cavity photon effec- 
tive mass. Therefore the exciton dispersion around the 
exciton-photon coupling can be neglected, and one can 
use w a (k) ~ Lo a - 4(Ji + J 2 ). 

The above Hamiltonian can be easily diagonalized to 
give 



(19) 



kr 



1) l^fel ~ 0' an d the upper 
branch becomes photonic, that is |^jl~| 2 ~ 0, \Yu~\ 2 « 1. 

In figure (4) we illustrate these results and plot the 
upper and lower polariton branches beside the cavity 
photon and exciton dispersions using the following typi- 
cal numbers. The exciton energy around small in-plane 
wave vector is taken to be hw a — 4(Ji + J2) = 2 eV. The 
distance between the cavity mirrors is L/m — 3100 A, 
which is chosen to give zero detuning between the exci- 
ton and the cavity photon dispersions at zero in-plane 
wave vector. For example for m = 3, we have iwl fim. 
The transition dipole is fj, = 2 eA, and the optical 
lattice constant is a = 2000 A. The exciton-photon 
coupling parameter, by using S = Na 2 , is given by 
l/l = yJ(u c (0)n 2 )/{2hLa 2 e ), with the value of |/| = 
4 x 10 11 Hz, where we neglected the k dependence for 
small in-plane wave vectors. At zero in-plane wave vector 
the Rabi splitting frequency of 8 x 10 11 Hz is clear. At 
large wave vectors the upper polariton branch tends to 
the cavity photon dispersion, and the lower branch tends 
to the exciton dispersion. 

In figure (5) we plot the exciton and photon weights 
in the lower polariton branch, |A~| 2 and |1"~| 2 . At zero 
in-plane wave vector the polariton is half exciton and half 
photon. For large in-plane wave vectors the polariton be- 
comes much more excitonic than photonic. In figure (6) 
we plot the exciton and photon weights in the upper po- 
lariton branch, |A+| 2 and \Y+\ 2 . At zero in-plane wave 
vector the polariton is half exciton and half photon. For 
large in-plane wave vectors the polariton becomes much 
more photonic than excitonic. 

2.0025 1 1 1 1 1 



which exhibits two diagonal polariton branches with dif- 
ferent dispersions 



(20) 



where = y/ 5\ + \fk\ 2 and we defined the exciton- 
photon detuning by 5k — (w c (fc) — u> a )/2. The splitting 
between the two polariton branches at wave vector k is 
2Afc and the splitting at the exciton-photon intersection 
point, where 5k = 0, is 2\fk\ which exactly corresponds 
to the vacuum Rabi splitting. 

The polaritons thus are coherent superpositions of ex- 
citons and photons with operators 



Ak± — X, Bu 



where the exciton and photon amplitudes are 



XI 



± 



Afc =F 5k 
2Ak 



fk 



y/2A k (A k T5k) 



(21) 



(22) 



At the exciton-photon intersection point the polaritons 
are half exciton and half photon, that is |^^| 2 = 1 2 = 
1/2. At large wave vectors the lower branch becomes 
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FIG. 4: The upper and the lower polariton branches vs. in- 
plane wave vector k. The dashed-dotted line is for the exciton 
dispersion, and the dashed-dotted parabola is for the cavity 
photon dispersion. 

Let us now discuss some dynamical consequences of 
this coupling and discuss the possibility to generate Bosc- 
Einstein Condensation (BEC) of polaritons analogous to 
recent setups in semiconductors [l8| for Wannier-Mott 
cavity polaritons. This process could be started by ex- 
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FIG. 5: The exciton and photon weights vs. in-plane wave 
vector k, in the lower polariton branch. 
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FIG. 6: The exciton and photon weights vs. in-plane wave 
vector k, in the upper polariton branch. 



citing lower branch polaritons at high energies. In turn 
the polaritons relax along the lower branch toward the 
minimum energy at k = 0. If one can achieve sufficient 
accumulation of polaritons at k — a phase transition 
into the BEC can be expected. The phase transition is 
stimulated by the nonlinear polariton-polariton interac- 
tions originated of the atom saturations. The relaxation 
towards k = is induced mainly by polariton-polariton 
interactions and to some extend by the emission of on- 
site vibrational quanta, which corresponds to localized 
atomic excitation to a higher Bloch band. As we ne- 
glect direct atomic tunneling such on-site excitations are 
dispersion-less, so that no bottle-neck effect exists in 
the present system. In semiconductor microcavities the 
main relaxation mechanism is due to the scattering of 
polaritons by acoustic-phonons and as the lower polari- 
ton branch slope around the exciton-photon intersection 



point is larger than the sound velocity in semiconductor 
crystals, the polaritons are stuck near this momentum in 
a so called bottle-neck effect Q. This effect is consid- 
ered as one of the main obstacles for achie ving a BEC of 
polaritons in semiconductor microcavities [18j |. Here we 
are centrally limited by the short life time of the atomic 
excited states and the cavity photons. The atomic decay 
time could in principle be enlarged within a cavity to get 
an effective polariton life time longer than the polariton 
relaxation. One alternative would be the use of long lived 
(metastablc) atomic states and very high finesse mirrors. 



V. POLARTION ANALYSIS VIA LINEAR 
CAVITY TRANSMISSION SPECTRA 

Up to now we considered an ideal cavity with perfect 
mirrors. But in order to observe the system eigenmodes, 
we need to couple the internal cavity modes to the exter- 
nal world. We do this in using the standard input out- 
put formalism or quasi-mode approach fl9| , which is well 
proved for high-Q cavities. In this approach the internal 
cavity field modes and the external free radiation field are 
quantized separately and weakly coupled via the cavity 
mirrors. This leads to the cavity field damping (photon 
loss) but also allows to consistently calculate the rela- 
tion between the cavity input and output fields. To in- 
clude atomic spontaneous emission as an energy loss term 
we add an effective exciton damping phenomenologically. 
This finally allows to calculate the transmission T, reflec- 
tion R, and absorption A spectra for a given incident field 
as shown in figure (3). The details of the calculations are 
presented in the appendix, which yield 



T( W ) 
R(u>) 
A{u) 



jirn I AM I 2 
|B( W )P ' 

l- ? 7[A(q;)-A*(a;)]+7 2 |A(a;)| 2 
z 7c/ [AH-A*H] 

I^MI 2 



(23) 



where we defined A(w) = J2 r 1^1 /(^ — 7 — (lu — 
7l)/2, |D(w)| 2 = |l+ry A(ui)\ 2 , and the other parameters 
are defined in the appendix. 

Now we present plots of the transmission, reflection, 
and absorption spectra. We consider the system given 
in the previous section and we plot the spectra for dif- 
ferent wave vectors. Here we choose the following num- 
bers for the damping rates. The exciton damping rate is 
T ex = 1.5 x 10 9 Hz, the upper and lower mirror damping 
rates are ju = 1l — 7.5 x 10 10 Hz. In figure (7) we 
plot the transmission, in figure (8) the reflection, and in 
figure (9) the absorption. The peaks of the transmission 
and the absorption, and the dips of the reflection corre- 
spond to the two polariton branches, which are the real 
eigenmodes of the system. From the plots we conclude 
that large absorption is obtained for polaritons which are 
much more excitonic than photonic. High transmission 
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peaks and deep reflection dips are obtained for polari- 
tons which are much more photonic than excitonic. At 
zero wave vector, as the polariton is half exciton and half 
photon, the upper and the lower spectra are identical. As 
the wave vector increase the upper branch becomes more 
photonic than excitonic, and the opposite for the lower 
one. Hence, the line-width of the upper branch spectrum 
becomes wider, which equals the cavity line-width, and 
the lower one becomes narrower, which equals the atom 
line-width. As the absorption is of the atomic excita- 
tions, at large wave vectors absorption only of the lower 
branch exists. 
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FIG. 9: The absorption spectrum for different in-plane wave 
vectors (the used parameters appear in the text). 



cavity polaritons. The atoms are taken to be two-level 
systems, where we chose two-levels close to resonance 
with a cavity mode, and their optical lattice potentials 
have equal light shifts. We limited ourselves to the Mott 
insulator phase with one atom per site. An excitation can 
transfer among the lattice atoms due to dipole-dipole in- 
teractions. The lattice symmetry allow us to present an 
excitation transfer by a wave that propagates in the lat- 
tice. These coherent states are Frenkel like-excitons, as 
in molecular crystals. We showed that excitons can be 
observed in optical lattices only if the atom line width is 
smaller than the exciton band width, which is achieved 
for optical lattices with a small lattice constant. After 
that, the system is fabricated between cavity mirrors. 
The excitons are coupled to the cavity photons by the 
electric dipole interaction. In the strong coupling regime 
the system eigenmodes are polaritons, which are coher- 
ent superpositions of excitons and photons. The polari- 
tons can be observed directly by linear optical spectra. 
The transmission, reflection, and absorption spectra for 
an incident external held show resonances at the polari- 
ton energies. The system presented here will have big 
applications in optoelectronic devices and quantum in- 
formation, and can open new directions in optical lattice 
systems, both theoretically and experimentally. Also the 
polaritons can be used as a strong observation tool for 
many physical effects in optical lattices. 



FIG. 8: The reflection spectrum for different in-plane wave 
vectors (the used parameters appear in the text). 
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APPENDIX A: LINEAR OPTICAL SPECTRA 
USING THE INPUT-OUTPUT FORMALISM 



and the polariton operator equation of motion is 



Here we calculate the linear optical spectra by adopting 
the input-output formalisms. The input and output fields 
are given by 

Hu = / dui k hw k b^(u k )bk(uk), 
k J 

Hl = ^ / dui k hw k cj c (w k )c k (wk), (Al) 
k J 

where ^(wk), ^k(^k) and c^(wk), Ck(^k) are the cre- 
ation and annihilation operators of an external mode 
with in-plane wave vector k and frequency Wk, above and 
below the cavity, respectively. The coupling between the 
external and internal fields at the cavity mirrors are given 
by 

k 3 

Vl = ^ / du k ihv(u k ) 4(cj k )a k - a^Ck(u! k ) 
k J 

(A2) 

where u(u>k) and i>(to>k) are the coupling parameters at 
the upper and the lower mirrors. The coupling is be- 
tween external and internal photons with the same in- 
plane wave vector. Here we consider light with a fixed 
polarization. For a fixed in-plane wave vector, we drop 
k, to get the total Hamiltonian 

h = ^2nn r a\a t 

duj hu b^(cj)b(uj) + J duj ftw c\uj)c{uj) 
+ f dwJ2 *M W ) [Y r * b\u)A r -Y r A\b{Lo)] 

J r 

+ f dujJ2 [Y r * J(u)A r - Y r Atc(cj)] , 

(A3) 

where, using the inverse transformation of Eqs. (f2~lj) , the 
cavity photon operator is written in terms of polariton 
operators a = '^2 r Y r * A r . The external radiation field 
operator equations of motion are 



dt 



b{uj) = -iw b(w) +u(uj)^2Y r * A r , 



d f 

— A r — —iQ, r A r — dui Y r \u(u>) b(u) + v(u) c(oj)] . 
dt J 

. , . (A5) 

The above system of equations are solved by using the 
input-output formalism [19(. The external field opera- 
tor equations are solved formally in using initial and fi- 
nal conditions, which are substituted back in the polari- 
ton operator equation and evaluated using a standard 
Markov approximation. For the upper mirror we define 
7u = 2itu 2 (uj) 1 and for the lower one 7^ = 2irv 2 (uj), 
which are related to the mirror reflectivity as jl oc 
(1 — Rl). Finally this gives: 



—A r = -iQ r A r + Y r (VTu" bi n + \flZ c in ) - jY r a, 

—A r = -ifl r A r - Y r (Vtu" b out + VlE c out ) + -yY r a, 

(A6) 



where 7 = (ju + 7l)/2, and we defined the input fields 
a,i n , Ci 7ll and the output fields a out , c ou t- The boundary 
conditions at the two mirrors are y/ju a = bi n + b out , and 
^/jl a = Ci rl + c ut ■ Exciton damping is included by using 
a complex polariton eigenfrequency £l r = Q r — iT r , where 
r r = r ex |X r | 2 , and where T ex is the effective exciton 
damping rate. 

By taking the Fourier transform we solve these equa- 
tions to get 



i(Q r — lu) A r 
i(fl r — lo) A r 
\/lu a 



Y T [y/jU b in + \FH ^n) ~ 7^ 

- Y ' r (y/lu bout + \/lL Cout) + lY r a, 

bin + bout , \/lL a = C in + Cout- (A7) 



For single side pumping, £i n — 0, we then can directly 
express the output fields as a function of the input field, 
with average photon numbers defined by n,i n — (b\ n bin), 
n out = (blutbout), and n L out = (c^„ t c ollt ). The transmis- 
sion coefficient T and the reflection coefficient R are de- 



fined by R(cu) = n^ ut /n in , and T(w) 



n out/ n i 



The 



j t c{uj) = -iuj c{uo)+v{uo)^Y r * A r . (A4) 



absorption, A, is calculated from the relation R(u>) 
T(uo) + A(u) = 1. The results are given in Eqs.(|23]), 
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